Weighing the spatial and temporal fluctuations of the dark universe 
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A generic prediction of the standard cosmology, based on general relativity (GR), dark matter and 
the cosmological constant (and more generally, smooth dark energy), is that, the two gravitational 
potentials describing the spatial and temporal scalar perturbations of the universe are equivalent. 
Modifications in GR or dark energy clustering in general violate this relation. Thus this ratio serves 
as a smoking gun of the dark universe. We propose a method to extract this ratio at various 
cosmological scales and redshifts from a set of measurements, in a model independent way. The 
ratio measured by future surveys has strong discriminating power for a variety of dark universe 
scenarios. 
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Introduction. — Predictions based on general relativ- 
ity (GR) plus the Standard Model of particle physics are 
at odds with a variety of independent astronomical ob- 
servations on galactic and cosmological scales, implying 
failures in particle physics or GR. There are various as- 
trophysical tools to probe this dark side of the physical 
universe (e.g. [H, Q)- Combining them allows us to break 
parameter degenaricies, reduce statistical errors and di- 
agnose possible systematics. 

These multiple probes are also crucial to detect smok- 
ing guns of new physics. For example, combining probes 
of the expansion history of the universe and probes of the 
large scale structure, the relation between the expansion 
rate and structure growth rate can be checked for signs of 
deviation from GR Indeed, one of the key questions 
in physics today is whether new particles/fields, such as 
dark matter and dark energy, or modifications to GR are 
needed to explain the observations. 

On large scales, two features of gravity can distinguish 
between a dark sector and modified gravity [1, [^, @, 0, [1] . 
One is the effective Newton's constant Goff, which spec- 
ifies the coupling between gravity and matter. In GR, 
Geff is equal to Newton's constant, but modified gravity 
models often predict deviations. The other is the rela- 
tion between the two gravitational potentials (p and ip. 
Here, the two potentials are defined in the Newtonian 
gauge through ds^ = (1 + '2tp)dt'^ — a^(l -f 24>)dy? where 
a{t) is the scale factor. The ratio r] = —(p/ilj weighs 
the relative ability of perturbations in matter-energy to 
distort the space-time.^ The standard cosmology, based 



on GR, dark matter and the cosmological constant (and 
more generally, smooth dark energy), predicts 77 = 1. 
Modifications from GR or emergence of intrinsic viscos- 
ity in dark energy fluid generally lead to 77 deviating from 
unity. Therefore, identifying observations, or sets of ob- 
servations, that will measure Goff and rj is of paramount 
importance 0, i, S, S 0, [HI • 

In we showed how to isolate the first key feature, 
feasibly testing the Poisson equation at ^ 1% accuracy 
level by combining weak lensing with galaxy redshift dis- 
tortion. In this paper, we will show that the same surveys 
allow us to directly measure 77, the second key feature, at 
cosmological scales. This can be done in a rather model 
independent manner. 

Models with rj ^ — 1. — Here we consider three models 
which produce deviations from the standard prediction 

77== 1 (04-7/;== 0). 

Perturbations in the Dvali-Gabadadze-Porrati (DGP) 
model have been carefully studied [l^, [3 • For a flat 
DGP model, 77 = [1 - 1/3/3dgp]/[1 + 1/3/3dgp], where 
/3dgp 1 - 2rc77(l + H/iH^) < and is the Hubble 
expansion rate. Here is the cross-over scale beyond 
which higher dimensional effects become important. In a 
flat model with matter density fc = l/^^o(l ~ ^m)- 
Since /3dgp < 0, ?/ > 1 in this model and the deviation 
from unity can be significant (Fig. [1]). 

Another modified gravity model (which aims to elimi- 
nate dark matter, not dark energy) is TeVeS [l^, a rel- 
ativistic version of MONd[i^. Besides the gravitational 
metric, TeVeS contains a scalar ps and a vector field. It 
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rj is the PPN parameter 7 (by forcing tli = —GM/r for point and sub-cluster scales are consistent with GR too [sl. [Tot. 
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has been shown [13, [Til that the TeVeS vector field can 
source the evolution of cosmological perturbations and 
compensate for the lack of dark matter in the model. 
To fit observations, the TeVeS parameter Kb should be 
small, in which case the vector perturbations a and E 
become large. These vector per turbations then drive ij 
to deviate from unity [l7, 18, T^ . 
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Here ( = {e~^'^^ — l)a. Since the background value (ps 
1 as imposed by nucleosynthesis bounds, the deviation of 
77 from unity is mainly driven by the vector perturbation 
(0 + "0 oc </>sa). A numerical evaluation of 77 is shown in 
Fig. [TJ For this figure we adopted a model with flf, = 
0.05, fl^ = 0.17, Oa = 0.78 and no dark matter. 

A final possibility is that gravity is still GR, but dark 
energy has non-negligible anisotropic stress a and causes 
inequality in two potentials through [2[ij 



(2) 



Although quintessence models predict a = 0, there arc 
some dark energy models that predict a and t] ^ 1 



2M,\2M- a specific example, we consider an extrinsic 
shear stress of the form rj = 1/(1+0;), with ui = a;oa"^(l — 
ilm)/f^Tn with ujQ constant, following Q. In general, ij 
varies not only with time, but also with scale. Richer 
physics encoded in the scale dependence of 77 would allow 
better discrimination between such dark energy model 
from other scenarios. 

The 7] estimator. — To measure 77, two independent 
measures of gravitational potentials are required. Both 
Vijj and V0 source the particle acceleration. However, 
the contribution from (j) is suppressed by a factor w^/c^, 
where v is the particle velocity. For this reason, non- 
relativistic particles such as galaxies only respond to ip. 
For the same reason, photons respond equally to both 
the potentials. Thus gravitational lensing measures the 
projected V^(0 — ip) along the line of sight. We propose 
an estimator consisting of the cross-correlation of each 
(the lensing field and the velocity field) with the galaxy 
distribution. 

The first cross-correlation is the lensing measure- 
ment with galaxy over-density in a narrow redshift 
bin[^. We can then obtain the cross-power spectrum 
Pv2(0-i/;)g(fc, z) between V^((/) — -ip) and the galaxy num- 
ber overdensity in the redshift bin associated with the 
galaxies. 

The second cross-correlation power spectrum Pgg can 
be obtained from the redshift distortions of the galaxy 
distribution in a spectroscopic survey 0, 0, IH, [23| . Here, 
6g = —V-Vg/H and Vg is the comoving peculiar velocity. 
We show below that this cross-spectrum is directly re- 
lated to Py^^/igj but first let us assume that this is so, that 
Pv'^ipg can be extracted from the 0-g cross-correlation. In 
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FIG. 1; Projected errors on 77 from SKA (half sky coverage is 
adopted). Solid line is the prediction of the standard ACDM, 
on which the error forecast is based. We also show the pre- 
dicted 77 for a flat DGP cosmology with Q,m = 0.2 (red, dotted 
line), the anisotropic shear stress model with loq = —0.3 (ma- 
genta, long dashed line) and TeVeS with Kb ~ 0.08 (blue, 
short dashed line). 



that case, the ratio of these two cross-spectra leads to an 
estimator for 77: 



P. 



v2(V'-?i)s 
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To see that 9 is related to Vi recall that on large 
scales gravity is the only force accelerating galaxies, so 
d{avP)/dt = — Vt/j, where = avg is the proper motion. 
Taking the divergence of this leads to 

^2^ ^ = _ {ln[a^HDg])' a^HV ■ Vg{a) . 

(4) 

Here, = d/da and Dg is the growth factor of 9g. The 
last relation holds in the linear regime where different 
modes decouple. We then have 



P^2^g = - {Ha^HDg]) a^H^Pi 



eg, 



(5) 



the desired relation. 

The proportionality factor relating the two cross- 
spectra in Eq. ([5]) requires knowledge of the expansion 
rate H{z) and the growth factor Dg[z). We assume 
that the former can be measured by other means; in- 
deed our goal is to distinguish dark sector models which 
produce identical expansion histories. No such assump- 
tion is needed for the growth factor, because the same 
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survey that measures Pgg will also measure Pgg, which is 
proportional to Dg and thus measurement of Pgg in mul- 
tiple redshift bins can be used to recover (In a^HDg) (see 
the appendix for details). We adopt the minimum vari- 
ance estimator to estimate errors in the reconstruction 
of Pgg and Pgg d, [31 • This reconstruction adopts no as- 
sumption on galaxy bias, so it is less affected by possible 
stochasticity or scale dependence in galaxy bias. 

Application of the 77 estimator in Eq. ([3]) relies on the 
condition of linear evolution such that Eq. ^ and there- 
fore ([5|) hold. For this reason, we restrict our discussion 
to the linear regime. This approach is robust against 
several uncertainties: (1) It docs not suffer uncertainty 
induced by the galaxy bias, whose effect cancels when 
taking the ratio in Eq. (2) It is not susceptible to pos- 
sible galaxy velocity bias, defined with respect to peculiar 
velocity of dark matter or dark energy, since we directly 
measure (hia^HDg) , instead of relying on a theory to 
calculate it. (3) It is applicable to general dark energy 
models and modified gravity models. It does not require 
dark energy to be smooth, nor gravity to be minimally 
coupled, nor scale- independent Dg. 

Forecast. — In order to measure 77 in this way, the lens- 
ing and redshift surveys must be sufficiently deep and 
wide. The proposed spectroscopic galaxy survey ADEPT 
or 21cm survey HSHS [2^, combined with a lensing sur- 
vey such as LSST, would be sufficient. Alternatively, 
SKA alone would be able to p rovide both suitable lens- 
ing, through cosmic shear [26| and cosmic magnification 
27|, and galaxy redshift measurements, as potentially 
would the Euclid^ mission. So we focus on SKA projec- 
tions. Dg can be measured by SKA at multiple bins of 
redshift and scale to impressive accuracy (Fig. [2]). We 
then infer {hi[a^HDg]) from the above measurements. 

Projections for the errors on 77 from SKA in a variety of 
{k, z) bins are shown in Fig.[TJ One example of the power 
of this measurement is in constraining the DGP model. 
The Eg measurement proposed in [ij can only marginally 
distinguish the f7„ = 0.2 fiat DGP model from ACDM. 
Fig.[T]shows, though, that these models have significantly 
different predictions for ?/; The TcVcS model adopted has 
been shown to produce a good fit of CMB and LSS data 
17| . However, with large deviation from 77 = 1, this 
model can be unambiguously distinguished from ACDM. 
Thus 77 and Eq are highly complementary to probe the 
dark universe'^; Modifications in gravity or dark energy 
viscosity often lead to stronger scale dependence in 77 than 
what is shown in Fig. [1] Our 77 estimator could have 
stronger discriminating power for these models. 
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^ http:/ /sci.esa.int/science-e/www/area/index.cfm?fareaid=102 
^ Errors in rj and in Eq are partly correlated. Future work should 

take this into account by fitting rj and Eq simultaneously, while 

marginalizing all other parameters. 



FIG. 2; Projected errors on De from SKA. Since SKA mea- 
sures Pee, which is proportional to Dg, this tells us directly 
on De, up to a normalization, which is of no physical impor- 
tance in rj measurement. The fiducial cosmology is the stan- 
dard ACDM cosmology (solid lines). The bin size is Ak = 
0.02/i/Mpc. For any realistic surveys, we have only limited k 
sampling. At radial direction, the available modes are fc, = 
2ttH {z)i/ Az, where i = 1, 2, • • • . This is the ultimate limiting 
factor of k sampling in wide field surveys such as SKA. Since 
we need to compare different redshifts to obtain ln(a^_De) , we 
require the k sampling at relevant redshift bins to be roughly 
identical. This requires Azi+i/H{zi+i) = Azi/H[zi). We 
choose {zi,Azi) = (0.4,0.2). Thus {z2,Az2) = (0.61,0.22), 
{z3,Az3) = (0.85, 0.26), [z4,A24] = [1.13,0.31], [z5,Azs] = 
[1.47,0.37], [z6,Az(i] = [1.88,0.45]. The maximum matter 
density power spectrum variance is max[A^(fc,z)] < 0.3 for 
the k and z bins shown, allowing us to neglect the effect of 
non-linearity for the moment. 



There is room to improve the 77 measurement. (1) The 
fractional error in Pgg reconstructed from the redshift 
distortion is about 15 times larger than the cosmic vari- 
ance limit [23 ]. Adopting the approximation of deter- 
ministic galaxy bias, the associated error will decrease 
by a factor of 3. However, to reach the cosmic variance 
limit, other velocity measurement techniques should be 
explored (e.g. |2g). (2) The forecast outlined above 
only uses Dg measurements from z < 2 galaxies and thus 
limited the accuracy of (\n.a?HDg) measurement. De- 
pending on the design and on the nature of 21cm emit- 
ting galaxies, SKA may allow measurements of Dg at 
higher redshifts. Furthermore, Dg at even higher red- 
shifts (z ^ 10) can be measured from redshift distortions 
of diffuse 21cm background. Improvement in the mea- 
surements of {Ina^HDg) and 77 at z < 2 would result 
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from the inclusion of such observations. (3) Furthermore, 
measurements of 77 at z > 2 can be made feasible by the 
inclusion of CMB lensing and 21cm background lensing. 

We have shown that future precision imaging surveys 
of weak gravitational lensing and spectroscopic surveys of 
galaxy redshift distortions provide highly complementary 
methods to probe the dark universe. In combination they 
allow us to isolate two key features of the dark universe, 
the effective Newton's constant Goff and r] = — from 
many astrophysical complexities, and distinguish com- 
peting scenarios of the dark universe robustly. 
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Appendix. — To infer {Ina^HDg) from Dg measured 
in limited redshift bins, a parametrization of Dg is re- 
quired. Since Dg evolve smoothly, {In a^HDg) should 
not be strongly dependent on the precise form of the 
parametrization. In this paper, we extend a widely used 
parameterization for Dg in standard gravity. For gravity 
models minimally coupled to matter, Dg — D a = JD, 
where / = d\nD/d\na and D is the linear density 
growth factor. One approximation adopted in the litera- 
ture is / ~ {nr,^a-^ I E^)'' (e.g. (H). Here, E = H/Hq is 
the normalized Hubble parameter. This approximation 
works well not only for CDM (7 = 5/9 for JIa + f7„j = 1 
and 7 = 0.6 for n\ = 0), but also for some modified 
gravity models such as DGP (7 = 2/3, We thus 

propose to fit a parameterization 



(6) 



Here, both Jl* and 7* are parameters to be fitted for 
each k bin. D and Dg are then obtained by the relation 
/ = d\nD/dlna and Dg = fD. 
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